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ABSTRACT 

Using atomic and molecular gas observations from the GASS and COLD GASS surveys and com- 
plementary optical/UV data from SDSS and GALEX, we investigate the nature of the variations in 
the molecular gas depletion time observed across the local massive galaxy population. The large and 
unbiased COLD GASS sample allows us for the first time to statistically assess the relative impor- 
tance of galaxy interactions, bar instabilities, morphologies and the presence of AGN in regulating star 
formation efficiency. We find that both the H2 mass fraction and depletion time vary as a function 
of the distance of a galaxy from the main sequence traced by star-forming galaxies in the SFR-M* 
plane. The longest gas depletion times are found in below-main sequence bulge-dominated galaxies 
(^*> 5 x 10 8 A1q kpc -2 , C > 2.6) that are either gas-poor (Mh 2 /M*<1.5%), or else on average less 
efficient by a factor of ~ 2 than disk-dominated galaxy at converting into stars any cold gas they 
may have. We find no link between the presence of AGN and these long depletion times. In the 
regime where galaxies are disc-dominated and gas-rich, the galaxies undergoing mergers or showing 
signs of morphological disruptions have the shortest molecular gas depletion times, while those host- 
ing strong stellar bars have only marginally higher global star formation efficiencies as compared to 
matched control samples. Our interpretation is that the molecular gas depletion time variations are 
caused by changes in the ratio between the gas mass traced by the CO(l-O) observations, and the 
gas mass in high density star-forming cores (as traced by observations of e.g. HCN(l-O)). While in- 
teractions, mergers and bar instabilities can locally increase pressure and raise the ratio of efficiently 
star-forming gas to CO-detected gas (therefore lowering the CO-based depletion time), massive bulges 
may prevent the formation of dense clumps by stabilizing gas disks against fragmentation, therefore 
producing the long depletion times. Building a sample representative of the local galaxy population 
with M*> 10 10 M Q , we derive a global Kennicutt-Schmidt star formation relation of slope 1.18 ± 0.24, 
and observe structure within the scatter around this relation, with galaxies having low (high) stellar 
mass surface densities lying systematically above (below) the mean relation, suggesting that £h 2 is 
not the only parameter driving the global star formation ability of a galaxy. 

Subject headings: Galaxies: evolution - Galaxies: ISM - Galaxies: star formation - ISM: molecules 



1. INTRODUCTION 

Star formation in galaxies proceeds at maxi- 
mum efficiency in systems undergoin g major mergers 
(|Genzel et al.l [20101: iDaddi et al.ll2010D . Given the well- 
established hierarchical structure formation paradigm, 
where larger structures form from the merging of smaller 
objects, it is tempting to infer that this "merging mode" 
of star formation contributes significantly to the build-up 
of stellar mass across the galaxy population, especially 
at high redshift. However, recent observations downplay 
the importance of violent events such as merger-induced 
starbursts in the glo bal star formation budge t of the uni- 
verse, even at z ~ 2 (|Rodighiero et alJl2Qll . 
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Instead, the existence of a tight relation between stellar 
mass and star formation rate, follo wed by star-forming 
galaxies up to at least z = 2 (e.g. iNoeske et al.l 120071 
Salim et all I2007t , IRodighiero et al.l I2010t INordon et al.l 



20101: lElbaz etld] 120111) 7 suggests that star formation is 



mainly regulated by secular processes with much longer 
duty cycles. To make further progress in understanding 
the histories of galaxies, we therefore need to investigate 
how processes such as minor mergers, bar formation and 
bulge growth influence the star formation process, both 
in terms of its onset / shutdown and of its efficiency. 

The first question to answer is thus, how do these 
mechanisms contribute to switching star formation 
on/off in galaxies? Several studies have recently tack- 
led this specific question, by looking for the observ- 
able that correlates best with quiescence. Both obser- 
vations and simulations suggest that quiescence is not 
strongly linked with stellar mass, but rather with struc- 
tural and morphological parameters such as Sersic index, 
stella r mass surface density, or central velocity disper- 
sion (IKauffmann et al.l I2003bt I Schiminovich et all 120071: 



More et al. 2011: Wuyts et al. 2011: Catinella et al.l 



20101 iCatinella et all l2012t iSaintonge et al.l l2011ab 
IKauffmann et al.ll2012D . While there is mounting evi- 
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dence for a link between star formation quenching and 
the presence of bulges, the direction of causality between 
them is a s of ye t not fully established. For example, 
iBell et all ()2011l ) suggest that while a bulge-dominated 
morphology is necessary for quenching to occur, it is not 
sufficient in itself. The exact underlying physical mech- 
anism^) responsible for the correlation between early- 
type morphologies and the absence of star formation has 
also yet to be unambiguously identified. We come back 
to this discussion in §4.31 

After having determined which galaxies are star- 
forming and which are quiescent, at least in a statis- 
tical way, the second question concerns the efficiency 
of star formation, across the galaxy population. The 
general picture is that a star-forming galaxy will con- 
vert a given amount of molecular gas into stars in a 
timescale of ~1 Gyr when star formation proceeds in 
self-gravitating Giant Molecular Clouds (GMCs), while 
in a major merger t his timescale is reduced by an or- 
der o f magnitude fe.g. lLerov et al . 2008; S olomon fe Sagd 
1988). It is however reasonable to expect that intermedi- 
ate systems exist, where the global (i.e. galaxy-wide av- 
eraged) star formation efficiency is found to vary within 
this f ramework (as suggested also by e.g. iMagnelli et al.l 
I2012D . Such situations could occur either in the onset of a 
major merger as the transition between the two modes of 
star formation occurs, or in systems subjected to milder 
processes where part of the star formation is occurring 
normally in GMCs, while at the same time more vio- 
lent star formation is taking place in other unvirialized 
regions. For ex ample, in the ov erlap region of the An- 
tennae galaxies. I Wei et all (|2012l ) observe a bimodal dis- 
tribution of molecular clouds, finding both a population 
of GMCs typical of nearby galaxies and a population 
of unus ually large, mas s ive, a nd efficiently star-forming 
clouds. iBoquien et al.l (|2011| ) also previously reported 
significantly different star formation efficiencies in dif- 
ferent regions of th e inter acting galaxy Arp 158, and 
iPapadopoulos et all (|2012j ) argue in general that accu- 
rate modeling of the interstellar medium (ISM) of lumi- 
nous infrared galaxies (LIRGs) requires consideration of 
both a central dense, efficiently star-forming phase and 
a more diffuse disc-like component where star formation 
proceeds at a more modest pace. 

Barred galaxies are examples of systems where the 
star formation process can be locally affected, as in- 
ward gas flows lead to increased central gas concentra- 
tions and star formation rates (e.g.|Sakamoto et al.lll999t 
ISheth et alJl2005tlMasters et al.ll2012fl . Interactions with 
nearby companions can also induce gas flows that cre- 
ate high levels of star formation in the central regions 
of ga laxies (e.g. iMihos fc HernquTstl 119961 : iKewlev et al.1 
2006). Conversely, some galaxies appear to resist star 
formation, even in the presence of cold gas. This in- 
cludes early-type, bulge-dominated systems, where the 
presence of large reserv oirs of atomic gas in the form 
of di sks is not rare (e.g. iGrossi et al.l l2009t iSerra et all 
1201 ll and references therein). Some of these early- 
type galaxies even have si gnificant, mostly compact, H? 
disks (e.g. iWiklind et all 119951: iKnapp fc Rupenl H996E 
lYoung et al.H201lD . suggested to be less efficient at form- 
ing stars than similar gas disks in normal star-fo rming 
galaxies (|Martig et al.l [20091: iCrocker et al.l I2012D . In- 
deed, several resolved studies point out that conditions 



in the ISM at the very centers of galaxies are different 
than in the disks, implying that the star formation effi- 
ciency is affected by the high st ellar densities and / or the 
presence of AGN dBoker et al.l[20TTt iMartm et al.ll2012b 
iDonovan Mever et al.ll2012l ). 

Here, we further investigate the link between structural 
parameters and star formation efficiency by looking di- 
rectly at the cold gas component of galaxies. Star forma- 
tion is a local process, and is thus optimally studied with 
observations at the highest resolutions possible. Our ap- 
proach here is however different: by working with inte- 
grated galaxy properties and gas measurements, which 
can be measured in large unbiased galaxy samples, we 
can track the processes that have an influence on the 
global star formation efficiency of galaxies, and there- 
fore the build-up of stellar mass at the present epoch. 
These results can then be linked to high resolution stud- 
ies of specific objects, to better understand the small 
scale physics responsible for the star formation efficiency 
variations. We use data from the COLD GASS sur- 
vey (CO Legacy Database for the Gafex-Arecibo-SDSS 
Survey: ISaintonge et al.l [201 lal ). which was specifically 
designed to conduct such statistical analyses. Com- 
pared with other samples which, w hile of comparable 
sizes, targeted specific populations (lYoung et al.l 119951 : 
ILisenfeld et all l20TTt lYoung et al.l [2011D, COLD GASS 
provides the full picture of the local galaxy population, 
over the mass range M*> 1O 1O M0. 

After describing the sample and the data in §21 we 
address in §3] t wo questions left open i n previous COLD 
GAS S papers (jSainton ge et al.ll2011bt iKauffmann et al.l 
120121) : (1) what is causing the molecular gas depletion 
time to vary across our sample? (2) what is the role 
of bulges in regulating and/or quenching star formation? 
As a sensor of the instantaneous star formation potential 
of galaxies, cold gas observations are a straightforward 
quantity that will help us answer these questions. A 
general interpretation of these results is presented in §3) 

All re st-frame a n d der ived quantities in this work as- 
sume a iChabrierl (|2003[ ) IMF, and a cosmology with 
H = 70km s" 1 Mpc" 1 , fi m = 0.3 and Q A = 0.7. 

2. SAMPLE AND DATA 

We have been conducting a full suite of surveys to 
characterize the stellar and gaseous contents of massive 
galaxies, in a sample selected to be representative of 
the entire galaxy population in the local Universe with 
M*> 10 10 M Q . The sample is selected in the redshift 
range 0.025 < z < 0.050 from the SDSS DR6 spectro- 
scopic survey, i n the area of sky t hat overlaps with the 
GALEX M IS (iMartin et al l [2005 1) and the ALFALFA 
HI survey ([Giovanelli et al.l 1200a lHavnes et~al~l 120 111 ). 
Full details regarding t he sample selection are g iven in 
iCatinella etall (fmof ): ISchiminovich et al.l (f2lnoh . The 
GALEX and SDSS DR7 photometry was reprocessed 
using a matched-aperture technique, in order to de- 
rive accurate colors, star formation rates, and stel- 
lar masses (see ICatinella et alJl20Tot iWang et all [2010: 
iSaintonge et al.l l2011al lbl). The star formation rates 
(SFRs) are obtained by fitting the photometry in the 
seven opt ical- UV bands to a s e t of m odel galaxies derived 
from the iBruzual fc Chariot! (|2003l ) code with a range 
of ages, star formation histories, metall icities and dust 
attenuations. For further details, see ISaintong e et al.l 
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(|20lTbl) . 

We also use in the analysis two structural parameters 
often derived from SDSS data to quantify the morphol- 
ogy of the galaxies. The first is the stellar mass surface 
density, defined as: 



where R^o.z is the radius in kpc encompassing 50% of the 
Petrosian flux in the SDSS z— band. The second param- 
eter, concentration index, is defined as C = -Rgo,r/^50.r- 
It has been shown to correlate well with bulge-to-total 
ratios as dete rmined by two-dimensi onal bulge/disk de- 
compositions ([Weinmann et al1l2009() . 

Out of the parent population of galaxies matching 
these selection criteria, a subset was ra ndomly selected 
for co ld gas observations. As detailed in Catincll a et al.l 
atomic gas measurements are obtained with the 
Arecibo telescope. The molecular gas masses come from 
observations of the CO(l-O) emission line with the IRAM 
30-m telescope. The sample used in this study includes 
the 222 galaxies that were part of the first COLD GASS 
data release, and the additional galaxies observed until 
the completion of the survey in 2011 November, for a to- 
tal of 365 galaxies. A full description of the instrumental 
setup, observi ng strategy and data red uction procedure 
i s pre sent ed in ISaintonge et al.l (|2011aD . We describe in 
5 ^2.21 and 12.31 the aspects of the process that differ from 
those used in previous COLD GASS studies. 

2.1. Sample definitions 
2.1.1. The full sample 

The COLD GASS sample is shown in Figure [TJ in the 
SFR-Af* plane. Since the sample is purely mass-selected, 
it includes both star-forming and passive galaxies. In 
order to quantify gas scaling relations across the entire 
galaxy population, the sample selected for gas observa- 
tions (both HI a nd CO) is constructed to have a flat 
di stribution ([Catinella et al.l I2010t ISaintonge et al.l 
l2011al ). Additionally, in order to fully explore the SFR- 
M* plane, we have targeted objects with high SFRs per 
unit mass. In this study, we refer to the 365 galaxies 
with CO measurements as the full sample (at the mo- 
ment, HI data are available for 320 of these galaxies). 
The full sample is meant to sample uniformly the entire 
SFR-Af* plane, and as such has an excess of high mass 
galaxies, and of galaxies with high specific star-formation 
rates compared to a purely volume-limited sample. 

2.1.2. The representative sample 

We also construct a representative sample, by select- 
ing a subset of objects from the full sample that matches 
the distribution of the unbiased, volume-limited parent 
sample in the SFR-Af* plane. We adopted a grid in this 
2-dimensional plane with cells of size 0.125 dex in M* 
and 0.2 dex in log(SFR), and culled galaxies from the 
full sample so that the relative fraction of the sample in 
each grid cell matched that in the parent sample. The 
distributions of other important quantities such as con- 
centration index and stellar mass surface densities, are 
also statistically equivalent between the parent sample 
and this representative subsample (Kolmogorov-Smirnov 
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Figure 1. Distribution of the sample in the SFR-Af* plane. The 
contours show the full unbiased SDSS parent sample or 12006 
galaxies, and the filled colored circles are galaxies with CO(1-0) 
measurements. The magenta symbols represent the galaxies in 
the main sample, the blue symbols the high SFR/M* galaxies. 
We will refer to all galaxies with a CO measurement as the full 
sample. We also define a representative sample, by picking a sub- 
set of the full sample that matches the distribution of the par- 
ent sample in this SFR-M* plane. An example of a representa- 
tive sample is shown by the open orange circle s. In this fi g ure w e 
also show the SFR-Af* relation determined by Peng et al. (120101 ') . 
log(SFR) = 0.9(0.1 log A/* - 1.0), and we define by the different 
lines six subsamples based on their distance from this star forma- 
tion main sequence. 

probability, KS pro b, greater than 0.8 that the C and /z* 
distributions are drawn from the same parent sample). 
We construct a set of 50 representative subsamples, an 
example of which is shown in Figure [TJ In the remainder 
of this paper, we explicitly mention whether the full or 
the representative sample is used. When the representa- 
tive sample is used, we are able to assess the robustness 
of the results by comparing the significance of the results 
across the set of 50 repetitions. 

2.1.3. The control samples 

In fj3J we investigate different mechanisms that may 
affect star formation efficiency. Because gas fractions 
and the de pletion time depend o n several global galaxy 
properties (ICatinella et all l2010t iCatinella et al.l l2012t 
ISaintonge et al.11201 laU bTlTit is not trivial to find a proper 
reference sample for any sub-population we wish to study 
(for example, galaxies with strong stellar bars or host- 
ing AGN). In these cases, we define a control sample 
of galaxies matched in some quantities, such that the 
effects of the mechanism under study on the depletion 
time can be assessed, everything else being equal. Un- 
less stated otherwise, we match the control samples on 
NUV— r color and stellar mass surface density, using 
the technique described above. The control sample there- 
fore shares the same distribution as the sample of interest 
in the the 2D plane defined by NUV— r and //*, the two 
quantities upon which atomic and molecular gas contents 
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depend the most dCatinella et al.l 120101 : ICatinella et al.l 
12011 ISaintonee et alJl20l"Tair For this reason, we must 
redefine a control sample each time we investigate a new 
sub-population of galaxies. 

2.2. Aperture corrections 

To correct the CO fluxes measured in the 22" beam 
of the IRAM telescope to the total flux, we derive an 
aperture correction for each galaxy. In the first COLD 
CASS data release, we used a small sample of nearby 
galaxies with resolved CO maps to derive a mean rela- 
tion for the aperture correction as a function of -D25, the 
SDSS g— band optical diameter. This method has the 
advantage of showing the range of aperture corrections 
that are possible at any given disk size -D25, given the 
real variation in radial molecular gas profiles within the 
reference sample. However, it does not take into account 
the specific properties of each galaxy, including its disk 
inclination and light profile. 

In this paper, we therefore use a new prescription to 
correct fluxes, simila r in spirit to the method used by 
ILisenfeld et al.l (|201 If) . For each galaxy with a CO mea- 
surement, we create a model galaxy having an exponen- 
tial molecular gas distribution with a "half-light" radius 
corresponding to r50.sF.fi, the radius encompassing 50% 
of the star formation as measured in the SDSS/ GALEX 
photometry. This choice is based on the observation that 
CO and SFR distributio ns trace each oth er well in nearby 
star-forming galaxies (jLerov et al.ll2009fl. even in t he HI- 
dominated outer disc regions (|Schruba et alllMTh . The 
model is then inclined as the real galaxy, and we simu- 
late an observation with the 22" IRAM beam. The ratio 
between this simulated observation and the total flux of 
the model galaxy is adopted as the aperture correction. 
For the more face-on systems, these values are perfectly 
consistent with the previous estimates. Only in the most 
inclined systems do we find aperture corrections on av- 
erage 15% smaller. 

The largest source of uncertainty in this process is 
the assumption that the gas profiles are exponential. 
Resolved maps of CO in nearby galaxies indicate that 
this assumption is generally valid: molecular gas pro- 
files are observed to be exponential, with the CO and 
SFR scale lengths equivalent, both o f which being pro- 
portional to the stellar scale lengths (lYoung et al.lll995l : 
iNishivama et al.l 120011 : iLerov et all 12008ft. albeit with 
some scatter especia lly at small radii (jRegan et al1l2001t 
iBigiel fc Blitd l2~012l ). At fixed disk inclination, we ob- 
serve differences of at most 15-20% in the values of 
Mh 2 obtained with aperture corrections derived from 
the method described above, and aperture corrections 
measured from the azimuthally-averaged observed SFR 
profile of each galaxy. Lisenfeld et al. (2011) estimate a 
30% uncertainty on their aperture corrections obtained 
from a similar disk modeling technique (see their Figure 
3). Based on our tests and this independent estimate, 
we adopt a general uncertainty of 20% on the aperture 
corrections. Accounting also for the uncertainties in the 
measurement of the CO line fluxes, in the flux calibra- 
tion at the telescope, and in the conversion factor aco 
(see ij2.3|) . we estimate a global error of ~ 0.3 dex on our 
values of M# 2 , with the error budget la rgely dominated 
by the systematic uncertainty on aco (Saintong e et all 
l20llal) . 



2.3. Molecular gas masses 

The other difference with the procedure described in 
iSaintonge et al.l (|2011aD is the choice of the CO-to-H 2 
conversion factor, aco- We have previously adopted a 
Galactic value of 3.2 M Q (K km s _1 pc 2 ) -1 for all the 
galaxies in our sample, which did not take into account 
the contribution of Helium. We now include this correc- 
tion, and therefore use aco= 4.35 Af Q (K km s _1 pc 2 ) 
for normal star-forming galaxies. Among the full sample 
used in this study, there are however galaxies undergoing 
major mergers. In such cases, it has been suggested that 
a value of 1.0 M^(K km s~ x pc 2 ) -1 is more appropriate 
jSdo mon et al.lll997l ). and we follow this convention. 

To avoid biasing the results of this study, it is prefer- 
able not to tie the choice between the Galactic and 
"merger" aco to a visual classification. Instead, we 
wish to use a quantitative measure of the properties of 
the interstellar medium to make this choice. While in- 
frared luminosity is directly connected to merging ac- 
tivity in the most extreme cases (ultralu minous infrared 
galaxies and sub-millimeter galaxies, e.g. Solomo n et al.1 
llQ^IVeilleux et al.l[2002T: lEngel et al.ll2010lV there is no 
one-to-one correlation between a high infrared luminos- 
ity {Lfir > 10 11 L Q ) and the need for a low value of 
acoi n or between morphological signs of merging and 
enhanced star formation ac tivity (e.g. iBushouse et al.l 
19881: lAlonso-Herrero etall 120061 : iDasvra et all 120061 : 
Cox et al.l 120081 ). However, the high densities found in 
mergers may lead to bot h elevated dust temper atures 
and lower values of aco ( Narayan an et all 1201 If ) . We 
therefore choose dust temperature, as parametrized by 
the Seo^m/ Sioo^m IRAS flux ratio, to guide our choice 
between Galactic and "merger" aco- Based on a sam- 
ple of galaxie s for which indepen dent measures of aco 
are available dBoselli et al.l I2002T ) . we choose a value of 
S&o^m/ SwOftm = 0.5 as the threshold to determine aco 
(Gracia-Carpio et al., in preparation). We therefore 
apply a conversion factor of aco= 1-0 M Q (K km s _1 
pc 2 )~ x to all galaxies with logLpiR/ 'L@ > 11.0 and 
S so fim / S 100 fim > 0.5, and the Galactic value of aco— 
4.35 M Q (K km s" 1 pc 2 )- 1 otherwise. Of the 365 galax- 
ies in the sample, 11 match this criteria and use the lower 
value of aco- We note that within our survey volume, 
IRAS was sensitive to galaxies with log LfirJLq > 10.7, 
therefore this convention to assign aco does not bias 
against part of the sample, even though only 15% of the 
COLD GASS galaxies have good quality flux measure- 
ments at both 60 and 100 fim. In fjB] we show that key 
results of this study are robust against this specific choice 
of a C o- 

3. RESULTS 

3.1. Depletion time variations across the SFR-M* plane 

We have already reported variations in both molecu- 
lar gas mass fraction (Jh 2 = Mh 2 / M*) and depletion 
time (t dep (E 2 )= M H JSFR) across the COLD GASS 
sample in the mass rang e 10.0 < log M^/M© < 11.5 
(jSaintonge et al.ll2011al fbl). These observations are sum- 
marized in Figure O We however improve here on our 
previous results by using stacking to calculate mean val- 
ues of £dep(H2) and fjj 2 , allowing us to include both CO 
detections and non-detections in the analysis. We com- 
pute these mean values in six bins of offset from the star- 
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Figure 2. Left: Mean depletion time (solid magenta line) and molecular gas fraction (dot-dashed black line) as a function of A(MS), 
the offset from the sequence formed by star- forming galaxies in the SFR-M* plane (Fig. [TJ. The mean value of td ep (H2) and fn 2 in 
each bin is obtained by stacking all CO spectra, therefore including both detections and non-detections, and the error bars obtained from 
bootstrapping arc lcr confidence intervals on the position of the mean. The location of the main sequence is indicated by the gray area. Both 
/jj 2 and the star formation efficiency (i.e. fd ep (H2) _1 ) increase in galaxies with high sSFR and decrease in passive galaxies. Comparing 
to the sSFR variation in the same bins (thick gray line), this figure shows that changes in sSFR across the SFR-Af* plane are not caused 
only by variations in gas contents, but also by changes in the star formation efficiency. Right: Mean concentration index (< C >, panel b) 
and stellar mass surface density (< log/i* >, panel c) in the same bins of A(MS). The black error bars show the lcr dispersion in each bin 
and the thick blue error bars the uncertainty on the position of the mean as determined by bootstrapping. 



formation main sequence, calculated as: 

A(MS) = log (SFR) -0.9(0.1 log M« - 1.0), (2) 



adopt ing the main sequence definition of iPeng et "ail 
( 2010). The different A (MS) bins are defined at fixed 
intervals, and illustrated in Figure [TJ Since the M*-SFR 
relation is almost linear , A (MS) is almost equivalent to 
variations in specific star formation rate (sSFR). 

To measure the mean depletion time of a given sub- 
sample, the stacking is done directly in "depletion time 
units" . That is, all multiplicative factors required to con- 
vert an observed CO line flux into a id ep (H2) value were 
applied to each spectrum prior to s tacking. This includes 
all the standard factors from the ISolomon et al.l (|1997T ) 
prescription (their equation 2), but also the appropriate 
value of aco (see H2.3[) . SFR -1 , and a weight to cor- 
rect for the flat stellar mass distributio n of our sample 
(the m ethodology for this is described inlSaintong e et all 
l2011al) . After stacking these scaled spectra, the mean de- 
pletion time is obtained by integrating over the resulting 
line profile, selecting by hand the spectral window over 
which the flux is integrated to reflect the effective line 
width in the stacked spectrum. This effective line width 
is seen to increase from ~ 300 to ~ 500 km s _1 going 
from low to high /i* or C, as these quantities correlate 
with M„, which in turn correlates with circular velocity. 

Similarly, we obtain mean values of fn 2 by stacking 
in "mass fraction units" (i.e. using a M* _1 scaling in- 
stead of SFR" 1 ). T his is the same technique used by 
iFabello et all (|2011al) to derive mean gas fractions from 
stacking, and further details can be found in that refer- 
ence. The only difference in methodology is that we do 



not weigh the spectra by their rms noise prior to stack- 
ing. Applying such weights would cause severe biases as 
the rms noise in the individual IRAM spectra is corre- 
lated with the underlying properties of the galaxies as 
a result of our survey strategy. Since we integrate un- 
til a galaxy is either detected in CO at S/N > 5 or a 
sensitivity to a gas mass fraction fn 2 =0.015 is reached, 
gas-rich galaxies and those at the smaller distances have 
spectra with higher rms than gas-poor and more distant 
objects. Applying a weight derived form this rms noise 
would therefore heavily bias against gas-rich objects. 

Figure [2]} shows that the mean value of both id ep (H2) 
and Jh 2 varies as a function of A (MS). Passive and 
weakly star-forming galaxies have low gas contents and 
long depletion times, while the most actively star- 
forming galaxies have higher gas contents and shorter 
depletion times. This is consistent for example with 
the conclusion of lCombes et al.l |l994) that the increased 
star formation rate in interacting galaxies is due to both 
increased molecular gas mass and increased star forma- 
tion efficiency. The COLD GASS sample shows over 
a much broader range of galaxy types that changes in 
the molecular gas mass fraction cannot alone account for 
the full amplitude of sSFR variations; compare the black 
dashed line and the solid gray line in Figure [2Ji. The 
mean star formation efficiency must also vary, and the 
stacking analysis indeed reveals that it increases by a 
factor of ~ 6 in the mean between passive galaxies in the 
red cloud and actively star-forming galaxies above the 
main sequence (see the magenta line in Fig. . 

What causes td ep (J^2) to vary between galaxies above, 
on and under the star formation main sequence? In 
€51 we argue that these changes are not caused by our 



6 



Saintonge et al. 



choice of the CO-to-H2 conversion factor. Studies of 
large optical datasets both at low and high redshift 
however indicate that other important physical param- 
eters vary smoothly as a func tion of distance from the 
star-formation seq uence (e.g. ISchiminovich et all I2007t 
IWuvts et alj|201ll ). and suggest that structural parame- 
ters quantifying bulge prominence are critical in setting 
the star formation properties of galaxies, or lack thereof. 
We confirm in Figure [2Jd,c that indeed our measures of 
"bulginess", the concentration index C and the stellar 
mass surface density do vary with A(MS), and in the 
next sections investigate the link between these struc- 
tural parameters and star formation efficiency variations. 

3.2. Depletion time variations with galaxy structural 
parameters 

The importance of structural parameters (or morphol- 
ogy) in setting the gas properties of galaxies has previ- 
ously been recognized. For example, there are thresh- 
olds in C and /i* above which the de tection rate of 
both HI and CO drops sign ificantly (jCatinella et al.l 
2012t ISaintonge et all 1201 laft . As also shown in 
Kauffmann et al.l (|2012l ). the C— /i* plane is particularly 
efficient at separating the quenched population from the 
gas-rich, star-forming galaxies. For this reason, we re- 
fer to the critical values of C = 2.6 and /z* = 10 87 M 
kpc~ 2 as "quenching thresholds". These critical values, 
determined from the gas properties of the galaxies, agree 
well with similar thresholds in stellar population or mor- 
phological indicators. For example, a concentration in- 
dex of 2.6 corresponds to the transition between the late- 
and early- type galaxy populations, as reported by stud- 
ies which compared this indicator to traditional visual 
morphologies, quantitative bulge/disk d ecompositions , 
and st ellar population indica t ors (e.g. |Strateya et all 
20011 IShTmasaku et all 120011: iNakamura et all 120031: 
Weinmann et all 120091: IKauffmann et alll2003bD . Simi- 
larly, a stellar ma ss surface density o f 3 x 10 8 M Q kpc -2 
was identified by IKauffma nn et all (|2003bl ) as marking 
the transition between galaxies dominated by young and 
old stellar populations. In comparison, stellar mass is 
a poor quantity t o separate the gas-rich from the gas- 
poor p opulations (|Catinella et al 1120101: ISaintonge et all 
I2011ah . 

In Figure[3l galaxies with cold gas mass fractions above 
the sensitivity limits of our surveys (f gas = 1.5%) are 
color-coded according to their £<jep(H2) (panel a) or their 
i(j e p(m) (panel b). This shows clearly that the galaxies 
with the longest values of tdepi^) and i^ep (HI) have high 
concentration indices and high stellar mass surface den- 
sities. It therefore appears that bulge-dominated, high 
stellar mass density galaxies either (1) have very little 
cold gas, cither atomic or molecular, or else (2) tend to 
be inefficient at converting into stars the molecular gas 
they do have. 

Care must be taken in interpreting this last observa- 
tion. Since at high ju* and C the detection rate of CO 
and HI is low, the galaxies we are detecting will tend to 
be the most gas-rich at fixed SFR, and therefore those 
having the longest depletion times. While it can be un- 
ambiguously concluded from Figure [3] that the longest 
values of tdep(J^2) and i<;ep(HI) are only found in bulge- 
dominated galaxies, this "Malmquist bias" means that 
we cannot conclude from Figure [3] alone that the mean 



depletion time is longer in bulge-dominated galaxies. To 
circumvent this issue, we can use stacking over represen- 
tative samples of galaxies with CO and HI spectra in 
order to include both detections and non-detections in 
the analysis. 

The mean values of tdepi^) and i<2ep(HI) are shown 
as a function of C and ju* in Figure [4] The stacking 
analysis confirms that both H2 and HI depletion times 
are longer in galaxies with high C and /i*, even when 
including the CO and HI non-detections. The specific 
behaviour of id e p(H2) and <d e p(HI) as a function of C 
and /i* is however different. The molecular gas deple- 
tion time increases steadily with C and /i» in the regime 
of disc-dominated galaxies located on or around the star 
formation main sequence (|A(MS)| < 0.4, see Fig. Hfc,d), 
and reaches a maximum value in the mean when reach- 
ing values of C and /i» characteristic of bulge-dominated 
galaxies. In contrast, td ep (HI) increases more slowly at 
C < 2.6 and log/i* < 8.7, but then suddenly reaches 
a value close to a Hubble time for the bulge-dominated 
galaxies, located well below the star formation main se- 
quence. An interpretation of these trends is presented 
below in i jOl 

3.3. Efficient star formation in dynamically active 
galaxies 

It is already established that in systems well above 
the main sequence of star-forming galaxies in the SFR- 
M* plane, tdepi^) is reduced by an order of magnitude 
compared to the normal star-forming population (e.g. 
Gao fc SolomoiJl2004HBothwell et aT]|2010t lGenzel et all 
20U[Daddi it al.l l2010). In the local universe, these are 
the ultraluminous infrared galaxies (ULIRGs), most of 
which (if not all) are majo r mergers ([Sanders fc Mirabel 
119961: iVeilleux et all 12002ft . For this reason, we spec- 
ulated in iSainton ge et al.l (|2011bl ) that milder dynami- 
cal effects, such as minor mergers and bar instabilities, 
could lead to more modest, yet measurable variations 
in the global idep(H2). We have now also established 
that tdep(H2) is longer in bulge-dominated galaxies than 
in star-forming disks, which may be interpreted in this 
same framework by an increased stability against star 
formation. In this section, we address the question of 
the impact of interactions and bars on star formation 
efficiency, and come back to the role of bulges in £14.31 

Whereas major mergers with mass ratios close to unity 
are very rare in the local Universe, the COLD CASS sam- 
ple includes a fair number of objects undergoing milder 
dynamical perturbations which may be affecting their 
star formation efficiency. The analysis here focuses on 
galaxies with bars, and on those displaying signs of in- 
teractions or morphological irregularities. We use the 
classification of Wang et al. (2012) to identify galaxies 
with strong bars, defined as having ellipticity etar > 0.5. 
Not only is the identification of such bars robust, these 
are the kinds of bars that are linked to nuclear starbursts 
(see Wang et al. 2012 for the full details). 

Two different approaches are then used to identify in- 
teracting or otherwise disturbed galaxies. First, we use 
a s et of quantitative m orphological parameters measured 
by iWang et al.l (1201 ID for t he GASS sample, following 
the technique of lLotz etaU (pOOl . Secondly, we per- 
formed a visual inspection of all COLD GASS galaxies, 
and classified them as either being normal, showing signs 
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Figure 3. Distribution of galaxies from the full sample in the /^*-C plane, color-coded by the depletion timescale of H2 (panel a) and of 
HI (panel b). In each case, non-detections are marked by open gray symbols, and the critical values of C and /i* above which most galaxies 
are non-detections are indicated by dotted lines. Above these values (i.e. in the top-right quadrant of each plot), galaxies either are very 
gas-poor, or tend to have long cold gas depletion times. 
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Figure 4. Molecular and atomic gas depletion times as a function of concentration index (a) and stellar mass surface density (b) in the 
representative sample. The mean values of i<iep(H2) and t,;,. p (0I) are obtained by stacking both detections and non-detections, and error 
bars are from bootstrapping (as in Fig. The histograms at the bottom of each panel shows the distribution of CO detections (solid gray 
histogram) and non-detections (dashed black histogram), and the shaded regions represent the ranges of C and fi t values where galaxies 
are mostly bulge-dominated and gas-poor. The detection rate as a function of fi* and C is similar for HI (Catinella ct al. 2013). In the 
bottom panels (c and d), the mean value of A (MS) is shown as a function of C and (i*. As in Figure [2] the black error bars are the la- 
dispersion in each bin and the thick blue error bars the uncertainty on the position of the mean as determined by bootstrapping. 



Saintonge et al. 



25 


representative 


sample 




25 


20 








20 


15 




Jl 




15 


10 








10 


5 


. _J 


hi 




5 







i_. . ■ 





full sample 













8.0 8.5 9.0 9.5 10.0 



8.0 8.5 9.0 9.5 10.0 




8.0 8.5 9.0 9.5 10.0 
l°9 t d . p (H 2 ) [yr] 



8.0 8.5 9.0 9.5 10.0 
l°9 t dsp (H 2 ) [yr] 



Figure 5. Molecular gas depletion time distribution for all COLD 
GASS galaxies with S/N C o > 5 and NUV-r< 5.0 (gray his- 
tograms). Top panels: galaxies undergoing mergers or showing 
signs of strong morphological disturbances (i.e. high asymmetries) 
are represented by the green histogram, while galaxies with strong 
bars are shown as the orange histogram. Bottom panels: depletion 
time distribution of the galaxies with p* < 5 X 10 s Mq kpc — 2 (blue) 
and /i*> 5 X 10 s Mq kpc — 2 (red). The left column shows results 
for the representative sample, while the right column displays the 
distributions in the full sample. 

of interaction and/or morphological disturbances, or fi- 
nally appearing to undergo a major merger. All galaxies 
identified as major mergers are IRAS-detected, with sig- 
nificant infrared luminosities (Lir > IO^Lq). In what 
follows we use the visual classification, but in gA]we show 
that not only there is a significant overlap between the 
subsamples identified by the two methods, the results of 
the analysis are not affected by the classification tech- 
nique. 

3.3.1. Interacting galaxies 

In Figure El the distribution of td ep (H.2) for the 
COLD GASS galaxies with reliable CO and SFR mea- 
surements is shown, for both the representative and 
full samples. In the top panels, the populations of 
interacting/morphologically-disturbed and barred galax- 
ies are highlighted. In the bottom panel of Figure [5] we 
show that a simple stellar mass surface density cut at 
log/it* = 8.7 is also efficient in separating the short- and 
long-depletion time populations in both the full and rep- 
resentative samples. The mean values of i<jep(H2) for the 
low and high /i* samples and the statistical significance of 
the difference between the distributions are summarized 
in TableQ] A KS test indicates that in both samples, low 
and high /i» systems have td ep (K2) distributions which 
have a probability < 0.001 of being equivalent. Since 
most of the morphologically disturbed galaxies have low 
stellar mass surface densities, to reliably assess the im- 
pact of interactions and bars on depletion time we define 
control samples of galaxies matched in fi* and NUV— r 
color but showing no signs of disturbance. 



As seen in Figure [5] and Table [TJ morphologically 
disturbed galaxies in the full sample have significantly 
shorter depletion times than objects in the control sam- 
ple (KS probability of 0.005 that th eir depletion times are 
drawn from the same distribution). iCombes et all (j!994f ) 
similarly observed the depletion time to be a factor of 2 
shorter in galaxy pairs. In gl we argue that the short 
tdep(R2) values in interacting/disturbed galaxies come 
from a change in the fraction of the gas that is found in 
the densest phase, due to the pressure increase. In the 
representative sample however, the depletion times of dis- 
turbed galaxies and their control are equivalent. This is 
consistent with previous findings that mergers cannot all 
be identicall y associated with episo des of active star for- 
mation (e.g. lDi Matteo et al.ll2007T) . either because some 
interactions may not drive gas inwards and compress it, 
or else because there is a delay between the occurrence 
of the signposts of morphological perturbations and star- 
burst events. 

3.3.2. Barred galaxies 

Figure [S] also suggests that barred galaxies have 
shorter-than-average depletion times. This is confirmed 
by the comparison with the control sample in both the 
full and representative samples, although the difference 
is small (a factor of 1.5), and the tdepi^) distributions 
only marginally different (KS pr . o &=0.25). Several stud- 
ies indicate that bars can drive gas inward and lead to 
locally enhanced nuc l ear star formation efficiency (e.g. 
ISakamoto et~aT1ll999t ISheth et all 120051) . In a detailed 
analysis of the barred population in the GASS parent 
sample, Wang et al. (2012) find a clear link between 
bars and increased central star formation in galaxies with 
M*> 3 x 10 10 M Q and ji»<3x 10 8 M Q kpc -2 , suggest- 
ing that bars are critical in the build-up of bulges, and 
possibly also in the quenching of star formation (see also 
iMasters et al.l 120111 ) . Our observations of only a mild 
increase of the global star formation efficiency in barred 
galaxies however suggests that bar instabilities do not 
significantly affect the star formation budget in the local 
Universe. Just as in the case of interacting galaxies, the 
week dependence of f,/ <7 ,(H2) on the presence of strong 
stellar bars may be due to the fact that we average over 
syste ms at different stag es of the gas inflow process. In- 
deed, [JogeeleF^I] (|2005l ) show that the depletion time of 
the molecular gas in barred galaxies can be both long or 
short, depending whether a galaxy is observed in early 
stages of the bar-driven inflow, or at a later stage when 
most of the molecular gas is in the circumnuclear region. 

3.3.3. The global star formation law 

To assess whether the tdepij^i) variations are caused 
by systematic changes in the global surface density of 
the molecular gas, the COLD GASS sample is shown 
in the Kennicutt-Schmidt (KS) plane in Figure |H1 Sur- 
face densities are computed assuming that st ar formation 
and m olecular gas are distributed similarly (jLerov et al.l 
2008), with a "half-light" radius, r 50 ^FR, measured for 
each galaxy as par t of our standard photometry pipeline 
(jWang et al.ll2010l) . Results are shown for both the full 
and representative samples. With the full sample one 
can assess the distribution of different galaxy types in 
the KS plane, while the representative sample allows the 
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Table 1 

Mean molecular gas depletion timescales a in selected COLD GASS sub-samples 







Bulges 






Bars 




Interactions 




Sample 


low fl* 


high )i t 




bars 


control 




interact control 


^-^prob 


Full 
Repr. 


0.96 ±0.44 
1.2 ±0.5 


1.9 ±0.9 
2.0 ±0.9 


< 0.001 

< 0.001 


0.88 ±0.35 
0.87 ±0.24 


1.3 ±0.8 
1.2 ±0.6 


0.26 
0.24 


0.63 ±0.41 1.1 ±0.6 
0.95 ±0.44 1.1 ±0.5 


0.005 
0.49 



a The mean depletion times are in units of 10 yr, and the errors quoted are the standard deviations around that mean 
value in any given subsample. 
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Figure 6. Kcnnicutt-Schmidt star formation relation from COLD GASS. Galaxies are color-coded by stellar mass surface density. The 
open symbols outline the position of "dynamically active" galaxies, and the dotted lines indicate constant depletion time values of 200 
Myr, 1 Gyr and 5 Gyr, from top to bottom. 



determination of an unbiased global star formation re- 
lation for local massive galaxies (M*> 10 10 M o ). This 
relation is: 



logS 



SFB, 



(0.86±0.21)(logS H2 ) - (2.94 ±0.26) 
(1.18 ±0.24) (log E H2 ) - (3.30 ±0.28), 



(3) 

where the errors given for the slope and intercept are 3cr 
formal fit errors, which approximates the total error in- 
cluding systematic effects and calibration uncertainties 
(|Genzel et al.l 12010). As different authors use different 
fitting strategies, we present for comparison purposes the 
results of an ordinary least squares regression of y on x 
(top equation), and of a bisector fit (bottom equation, 
and our preferred method). In both cases, the scatter 
of the residuals around the relation is <jrs = 0.28, and 
the slope is consistent with linear. These values are in 
agreement with oth er derivations of the star formation 
law at the kpc-scale (TLerov et al.l|2008HBlanc et alj|2009t 
iBigiel et al.l 1201 it iRahman et al.ll2012h and at high red- 
shift (|Genzel et al.ll2010UDaddi et alJl2010f >. 

The £<jep (Eb) variations observed across the COLD 
GASS sample do not come from a star formation law 
with significantly steeper slope than in other recent work, 
nor from systematic trends between X# 2 and e.g. jit* or 
interaction state. Rather, the idep(H2) variations mani- 
fest themselves as structure within the scatter around a 
^linear KS relation, where specific populations, for ex- 
ample galaxies with high/low /i* are preferentially lo- 
cated below/above the mean KS relation (Fig, [5]). As 
expected from previous studies, the galaxies undergo- 
ing major mergers have short depletion times and mark 



the upper envelope of the distribution in the KS plane 
(Fig. E}. Mildly interacting/disturbed galaxies, having 
low stellar mass surface densities, then fall systematically 
above the mean KS relation. 

Additionally, galaxies with the highest stellar mass sur- 
face densities and concentration indices (i.e. the bulge- 
dominated galaxies) are found systematically below the 
mean KS relation, with long global molecular gas deple- 
tion times ranging between 1 and 5 Gyr. A similar obser- 
vation was made bv lCrocker et al.l (|2012l ) for a small sam- 
ple of nearby early-type galaxies with detailed observa- 
tions of multiple molecular gas lines. In Figure UK, there 
are a few galaxies deviating from the rest of the popula- 
tion, having high stellar mass surface densities yet short 
depletion times. These are however all merging systems, 
which may be in a transition phase and will ultimately 
lie below the mean KS relati on after the starburst p hase, 
as suggested by simulations (|Bournaud et al.lf201 lh . We 
come back to this discussion in §4.31 

The main conclusions to be drawn from Figure [6] 
are that (1) in the representative sample, the "merger 
branch" of the KS relation mostly disappears, confirming 
that very efficiently star-forming major mergers are rare 
events that do not contribute significantly to the star 
formation budget in the local Universe, (2) the global 
surface density of H2 (estimated form our CO (1-0) line 
fluxes) is not the only parameter driving the star forma- 
tion efficiency of galaxies, as we observe structure within 
the scatter about the mean KS relation, and (3) the 
bulge-dominated galaxies, which have the longest molec- 
ular gas depletion times, are inefficient star-formers de- 
spite have Tih 2 values as high as spiral galaxies (see also 
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lYoung et al. l l200a IShapfro et al.ll20Toh . 

3.4. The role of AGN 

As explained in fJTJ most studies suggest that quench- 
ing is not strongly linked with stellar mass, but rather 
with a measure of "bulginess" . Given the well- 
established correl ation between supermassive black holes 
and b ulges (e.g. lHaring fc Rbd 120041: iMagorrian et al.l 
[1991 IMerloni et al.l 120101: IBeifiori et al.l 120121) . AGN 
feedback as an important quenchi ng mechanism is ap- 
pealing (e.g. ISomerville et al.ll2008l ). as it would explain 
why few galaxies without central mass concentration 
are quenched, while a very large frac tion of the bulge- 
dominated systems are (e.g. lBellll2008D . 

In this light, there has been work done to investigate 
the link between the presence of AGN and the gas content 
of galaxies, which sh ould be most dire ctly affected by the 
feedbac k. However, IHo et al.l (|2008 ) and iFabello et al.l 
(|2011bf) report no difference in the atomic gas properties 
of active and inactive galaxies, in large representative 
nearby samples. This would seem to argue against the 
AGN feedback scenario presented above, although both 
studies point out that molecular gas, which is generally 
significantly more centrally-concentrated than HI, may 
be a more sensitive tracer. Here, we use our IRAM CO 
observations to investigate whether the presence of an 
active nucleus influences the molecular gas contents of 
galaxies, in order to address the question of why bulge- 
dominated galaxies have on average longer td ep (H2 ) . As 
detailed below, we find that while AGN may affect the 
amount of molecular gas in their host galaxies, the star 
formation rate is also reduced, leaving the depletion time 
unchanged (or if anything reduced) as compared with a 
properly-matched control sample. 

We iden tify as active the COLD G ASS galaxies located 
above the Kauffma nn et al.l (|2003a[) line in the BPT di- 
agram (the line rat io diagnostic plot first presented by 
iBaldwin et al.|[l98l"l) . Of these, 6% are Seyferts, the rest 
being either LINERs or composite systems where part 
of the line emission is from star formation. In the refer- 
ence sample of "non-AGN" we include the star forming 
systems identified from the BPT diagram and the in- 
active galaxies which have no line emission (S/N < 3 
in Ha, [Nil], H/3 and [OIII]). To compare fairly their 
molecular gas properties, we match subsamples of AGN 
hosts and reference galaxies, requiring that the two sub- 
samples have equivalent distribution in /j,* and NUV— r 
color. The results presented below are unchanged if in- 
stead we control for star formation activity in the central 
region (i.e. using the fiber D„(4000) instead of the global 
NUV— r color), or for another structural parameters such 
as or C. 

As a proxy for the Eddington ratio, we use the ra- 
tio L[OIII] /ct 4 . The luminosity in t he [QUI] line whic h 
we correct for extinction following W ild et all ([2007), 
can be used as a proxy for the black hole accretion 
rate (fKauffma nn et af1l2003a|) . while the black hole mass 
is proportional to a 4 (|Tremaine et alj 12002). We use 
for a the stellar velocity dispersion measured in the 
SPSS fiber, correct ed for aperture effects as detailed in 
IGraves et al.l ([20091 ) . Each subsample (AGN and con- 
trol) is divided in three bins of _L[OIII]/(T 4 and galaxies 
are stacked in order to obtain the mean values of fu 2 



and idep(H2). The results are shown in Figure [7] At all 
values of the Eddington ratio, the mean molecular gas 
mass fraction is lower in the AGN subsample, but the 
mean depletion times are consistent or even shorter than 
i n the control. 

IMaiolino et al.l (|1997l ) conducted a similar study, com- 
paring the molecular gas properties of Seyfert and field 
galaxies. They report no strong difference in the global 
H2 contents, but an increased star formation efficiency 
in Seyfert 2 hosts. Our very tentative observation of 
a shorter tdepi^) in AGN hosts may be the extension 
of what these authors are finding in the more powerful 
Seyferts, although a direct test of this is not possible due 
to the lack of overlap between the types of AGN probed 
by our studies and the different methodology in assigning 
control samples. 

The most important conclusion of this section in the 
context of this work is that the high incidence of AGN in 
bulge-dominated galaxies cannot explain why tdepi^) is 
longer at high concentration index and There does 
however appear to be a link between AGN and fn 2 , but 
we caution that this observation pertains only to the type 
of AGN which are in our sample, typically LINERs ac- 
creting at a rate of at most ~ Q.QlLEdd- 

4. DISCUSSION 

4.1. Summary of observational results 

The key observations presented in this work can be 
summarized as follows: 

1. Both the molecular gas mass fraction (/h 2 ) an d 
depletion time (td ep (H 2 )) vary across the SFR-Af* 
plane. Galaxies above the star-forming sequence 
have both higher fn 2 and shorter trf ep (H2). Con- 
versely, explaining the low specific star formation 
rates of galaxies below the star-forming sequence 
requires both fn 2 to decrease and tdepi^) to in- 
crease. 

2. The bulge-dominated galaxies in the sample, iden- 
tified by their large concentration indices or their 
high stellar mass surface densities (C > 2.6, //*> 
10 8 7 M Q kpc -2 ) either have no cold gas down to 
the limit of our survey, or else are on average in- 
efficient at converting into stars the molecular gas 
they do have. Despite the fact that molecular gas 
in early-type galaxies has high surface densities 
just like in normal star-forming discs, we find their 
mean td e p(R2) to be twice as long. The atomic gas 
depletion time is also longer in bulge-dominated 
galaxies. 

3. At fixed /J>* and color, galaxies hosting AGN have 
lower molecular gas mass fractions than their con- 
trol sample, at all values of the Eddington ratio 
probed by our sample. There is tentative evidence 
that the depletion timescale of the molecular gas 
is either unchanged or lowered in the presence of 
an AGN, but not increased, indicating that active 
nuclei are not the main agent responsible for the 
long depletion times in bulge-dominated galaxies. 

4. Among the gas-rich, disc-dominated population, 
galaxies undergoing merging events or having dis- 
turbed morphologies have the shortest depletion 
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Figure 7. H2 mass fraction and depletion time as a function of the Eddington ratio, parametrized as L[OIII]/cr 4 . Individual CO detections 
and non-detections are shown, as well as mean values in bins of L[OIII]/<t 4 obtained by stacking (colored symbols, red: AGN, blue: control 
sample). Galaxies in the AGN subsample (filled circles) have lower molecular gas mass fractions than galaxies in the control sample 
(squares), but their depiction times are equivalent if not shorter. 



times, on average a factor of 2 shorter compared 
to matched control samples. Not all mergers or in- 
teractions are however associated with episodes of 
particularly efficient star formation. 

5. Global molecular gas depletion times are only 
mildly shorter in galaxies hosting strong stellar 
bars, everything else being equal. Even though re- 
solved gas studies show that bars can lead to en- 
hanced nuclear star formation efficiency the effect 
is not strong enough to significantly affect the global 
depletion time. 

6. We derive the star formation law for an unbiased 
sample of galaxies, representative of the local pop- 
ulation with M*> 10 10 M Q . This relation has a 
slope of 1.18 ± 0.24, consistent with previous stud- 
ies; the tdepi^) variations manifest themselves as 
structure within the scatter around this best-fit re- 
lation, with bulge-dominated galaxies falling sys- 
tematically below the mean relation. 

We now discuss these results, addressing the two main 
questions posed in the Introduction, and comparing them 
to previously published work. 

4.2. What is causing the global star formation efficiency 
to vary across the local massive galaxy population ? 

Our results confirm previous observations that the 
transition from star forming to quiescent galaxies is cor- 
related with the presence of a bulge, as quantified cither 
by stellar mass surface density, Sersic ind ex, velocity dis- 
persi o n, or concentration inde x (see also [C atinell a et al.l 
[20101: iSaintonge et~aH l2011at iKauffmann et al.l 12012). 
While there is still some debate as to which of these 
measures of "bulginess" correlates be st with quiescence 
(for different view s on this, see e.g . iFranx et al.l l2008t 
IWuvts et al.1 l2lHl IBell et al.l 12011 IWake et al.1 I2012T) . 
various studies generally agree that stellar mass is a 
poorer diagnostic to separate star forming from quiescent 
galaxies. Our results bring additional weight to these 



previous studies, since they are based on cold gas mea- 
surements, a direct tracer of the instantaneous potential 
of a galaxy to form stars, rather than on indicators such 
as color, which are averages over the different star for- 
mation history of each galaxy. 

Structural parameters such as /i* and C do not only 
correlate with quiescence, but also with star formation 
efficiency. These physical quantities are therefore at the 
core of this question. Only in galaxies with lower /1* do 
we see the dynamical processes that can significantly in- 
crease star formation efficiency. At the other extreme, 
high central concentrations may lead to reduced star for- 
mation efficiency or complete quenching of the star for- 
mation. 

The single quantity that may link these observa- 
tions is the dense molecular gas fraction. We specu- 
late that while td e p(H-2) (measured from CO(l-O) line 
fluxes) varies across our sample, were we to measure 
a depletion time for the dense gas, it would show sig- 
nificantly less scatter. A commonly used dense gas 
tracer is the HCN(l-O) line, with observations showing 
that the HCN(1-0)/CO(1-0) ratio is i ncreased in effi- 
cient ly star-forming galaxie s like LIRGs dGao fc Solomonl 
2004t Uuneau et al. I2009t iGracia-Carpio et al.l 120081: 
Garcia-Burillo et al.l I2012D. and possib ly reduced in 
early- type galaxies ([Crocker et al.l 120121 ) . We caution 
however that although the HCN/CO ratio appears to be 
a signpost for the star f orming dense gas fraction, it is 
not a perfect tracer (e.g. iGracia- Carpio et al.l [2006). In- 
deed, HCN is a minor species subject to the abundance 
of Nitrogen and the details of the chemical network (e.g. 
Sternberg et~aT1 Il994t lUsero et al.l 120041; IMeiierink et al.l 
20071) . and its (1-0) transition is sensitive to the specific 
physical excitation condition s with mid-infrare d pump- 
ing possibly playing a role (|Aalto et al.l 119951 ). Other 
possibilities to trace the dense star-forming gas would 
be the HCO + (1-0) transition, or higher rotational tran- 
sitions of e.g. CO, HCN or HCO+. 

Th eoretical and numerical work al so supports this pic- 
ture. IKrumholz fc Thompson! ([2007D find that a line with 
a high critical density like HCN(l-O) traces well the ef- 
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ficiently star forming gas between galaxies with differ- 
ent global properties. On the other hand, they find 
that a lower density tracer like CO(l-O) tracks regions 
with various star formation properties from galaxy to 
galaxy. With high reso lution hydrodynamic simulations, 
iBournaud et al.l (|2011h find the probability distribution 
function (PDF) of the local gas density in mergers to 
be broader and shifted to higher densities (or perhaps to 
have a different functional form) , compared to the mostly 
log- normal PDFs found in normal star-forming galax- 
ies (jWada fc Normanl 12001). They conclude that these 
variations in the dense gas fractions at different stages 
of a merger cause galaxies to go from having "normal" 
i(2ep(H2) before the interaction begins, to having short 
tdep{)^2) at the peak of the merging activity (resulting 
in the merger sequence in the Kennicutt-Schmidt plane). 
The end product of these mergers, early-type systems 
with some amount of dense molecular gas left, are some- 
what inefficient at forming stars and have longer than 
average td ep {)^2)- This scenario explains well the loca- 
tion of the various galaxy types on the COLD GASS KS 
relation (Fig. [6]). 

4.3. What is the role of bulges in regulating/quenching 
star formation? 

In early simulation work, IQstriker fc Peebles] (|1973T ) 
pointed out that galaxies with significant spheroidal com- 
ponents are more stable against fragmentation than pure 
disks. Since fragmentation introduced by gravitational 
instabilities is a key element of star formation, the nat- 
ural conclusion may be that bulges can therefore re- 
duce star formation efficiency. This can be understoo d 
in terms of the stability parameter, Q (|Toomrd 11964( 1. 
which in the case of gas+stars sys tems can be extended 
to include the two components (e.g. I Jog fc SolomoriT984l ; 
lElmegreenl 119951 : [Romeo fc Wiegertll2011f), as confirmed 
by observations (IWong fc Blitz] 120021: iLerov et "ail 120081: 
IHitschfeld etaLl |2009J). Specifically, a gas disk, which 
would be otherwise unstable, can be made stable by a 
significant stellar bulge. Under the framework presented 
above, long depletion times in bulge-dominated galaxies 
would therefore result from a gas density PDF shifted to 
lower densities than in normal star-forming galaxies. We 
note that the effect of a stellar component on the stabil- 
ity of a gas disk could however be reversed in galaxies 
with lower stellar masses than those in the COLD GASS 
sample, when stars are distributed in a disk without a sig- 
nificant central concentration. In that case, higher stellar 
mass surface densities could translate into a sho r ter de - 
pletion time, as argued for example bv lShi et al.l (|2011[) . 

As seen in Figure [6l bulge-dominated galaxies fall 
systematically below the mean molecular Kennicutt- 
Schmidt relation, indicating that at fixed gas surface 
density, they are less efficient at forming stars. The ob- 
servation is in line with the scenario presented above, in 
which massive stellar bulges can stabilize gas disks, there- 
fore reducing the star formation efficiency. Although 
this mechanism is not sufficient in itself to make star- 
forming galaxies "red and dead", it may be efficient at 
keeping already-quenched galaxies on the red sequence, 
by preventing further star form ation even in the e vent 
of late-time cold gas accretion (Ma rtig et al.l l2009h . In 
massive galaxies, radio-mode AGN feedback may be the 
most important agent in keeping galaxies on the red se- 



quence (e.g. lGabor et al.ll2011l:lvan de Voort et al.ll2011ft . 
but the contribution of the bulge in the reduction of the 
star formation efficiency could however be an important 
effect in lower mass galaxies where AGN are not preva- 
lent. 

Interestingly, we also find the longest values of i<2 ep (ffl) 
in galaxies above the "quenching thresholds" (C > 2.6 
and 10 8 7 M Q kpc -2 , see Fig. EJ, although in this 
case the reason may be different, and the different be- 
havior of tdep{)^2) and i(jep(HI) as a function of [i* hints 
at changes in the bottleneck of the star formation pro- 
cess across the local galaxy population. In high density 
bulge-dominated galaxies, the atomic gas depletion time 
increases significantly. This suggests that these galaxies 
have large reservoirs of HI gas, presumably located in the 
outer disk, that do not participate in the star formation 
process. In this regime, the bottleneck for star formation 
is the transport of this atomic gas to regions of higher 
density where the atomic-to-molecular conversion (and 
hence star formation) takes place efficiently. At stellar 
mass surface densities characteristic of disc-dominated 
galaxies, td ep (HI) has a mean value of ~ 3 Gyr with larg e 
galaxy-to-galaxy variations (|Schiminovich et"ail 120101) . 
but increases slowly with /z» and C similarly to td ep (H-2)- 
On average, there is therefore no accumulation of atomic 
gas, the conversion to CO(l-0)-detectable molecular gas 
proceeding unimpeded. The ^^(H^) variations then 
come about as galaxies have different fractions of their 
molecular gas in efficiently star-forming dense cores. 
This interpretation is supported by our observation that 
even among galaxies with low stellar mass surface den- 
sities, the ones undergoing merging, or to some extent 
those hosting a strong stellar bar, have the shortest CO- 
based ^ ep (H 2 ). 

Since bulge-dominated galaxies preferentially host 
AGN, it has been suggested that the correlation between 
quenching and the presence of a bulge may be a conse- 
quence of the feedback from the AGN. If this were the 
case, there should be differences in the gas contents of 
galaxies with and without AGN, everything else being 
equal. Compared to a control sample matched in /z* and 
color (either global NUV— r or central D„(4000)), galax- 
ies hosting an AGN with — 4 < log L \oui] /°~ A < — 2 have 
lower molecular gas mass fractions than inactive galax- 
ies. This observation indicates a possible connection be- 
tween AGN activity and the quenching of star forma- 
tion, although our sample does not include the stronger 
AGN which are likely to contribute most to this pro- 
cess. Galaxies with and without AGN in our matched 
subsamples do not have statistically different mean val- 
ues of idep(H2), however, suggesting that the AGN with 
modest accretion rates are not responsible for the long 
depletion times in bulge-dominated galaxies. There is 
very marginal evidence that td ep (^2) is even reduced in 
AGN hosts, although the difference is not statistically 
significant in our sample which incl udes mostly LINERs . 
The non-result is in agreement with Uuneau et al.l (|2009f l 
that report no difference in the HCN/CO ratio between 
galaxies with and without AGN, which in the framework 
of this discussion could explain why we do not observe 
variations in the CO-based values of td ep (H2) for AGN 
hosting galaxies. The disconnect between nuclear prop- 
erties and outer disk star formation explains why there is 
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no correlation between AGN activity an d the generally 
extended HI contents of the host galaxies (iHo et al.ll2008t 
iFabello et a!1l2011bl : lDTamond-Stanic fc Riekdl2012D . 

4.4. Comparison with results at the kiloparsec scale 

One of the main results of this work is that global struc- 
tural parameters, in this case concentration index and 
stellar mass surface density, are critical in establishing 
both the bimodality between gas-rich star-forming and 
gas-poor quiescent galaxies, and in setting the efficiency 
of star formation. Star formation however is not purely 
a global process, as it depends strongly on small-scale 
physical conditions, which our integrated measurements 
cannot track. Significant work has however gone into 
measuring gas and star formation in nearby galaxies at 
the (sub-)kiloparsec scale, with results, if taken at face 
value, that may appear at odds w ith the conclusions of 
this study (e.g. iBigiel et al.ll2011l who present evidence 
for a constant depletion time in nearby galaxies). 

The relation we observe between i<j ep (H2) and the 
global stellar mass surface de nsity may not appl y at small 
scales, although Figure 15 of lLerov et al.l (|2008l ) indicates 
that the local mean depletion time in kpc-sized regions 
of nearby spirals increases by a factor of two over one 
de cade in (a simila r effect can be seen in Figure 2 
of iRahman et alj 12012). The global values of /i* may 
not track the local conditions in individual sites of star 
formation, but rather the probability of incidence of var- 
ious dynamical processes which can influence whether 
star formation proceeds "normally" in virialized GMCs 
or rather in dense super star-forming complexes. 

In a study of th e resolved star formation law, 
ISchruba et all (l20TTh argue that star formation effi- 
ciency does not vary significantly within individual galax- 
ies, with however significant galaxy-to-galaxy variations. 
Even though this is consistent with the observations pre- 
sented here, the interpretat i on is fundamentally differ- 
ent. Indeed, ISchruba et al.l (|201lD ascribe the idepC^) 
variations to a metallicity dependence of the conversion 
factor aco, while our analysis suggests that they are 
caused by actual changes in star formation efficiency. 
As detailed in the appendix, we do not find any evi- 
dence for metallicity or extinction to be the main driver 
of the observed tdepi^) variations. In the stellar mass 
range of the COLD GASS sample (> 10 10 M W ), t he mass- 
metallicity relation is flat (jTremonti et al.l [20041 ) . imply- 
ing that quantities such as /i* and C do not correlate 
with metallicity. However, toward the low mass end 
of our sample, and even more so at M*< 10 10 M Q , the 
M# 2 /SFR ratio may very well be modulated simultane- 
ously by variations in aco and actual changes in star 
formation efficiency. Even after accounting for metal- 
licity effects on aco, high star formation efficiencies 
have indeed been reported for several low mass, low 
/it* local galaxies where high resolution measurements 
of both g as and star formation are possible, for exam- 
ple M33 (IGardan et al.l 120071: IGratier et al1l20Tol) IC10 
(jLerovet al.l 120061 ). and NGC6822 (|Gratier et al.l f2(TToh 
(see ho wever the counter a rgument in the case of the 
SMC in lBolatto et al.l 120111 ). Disentangling the two ef- 
fects will require both large statistical samples of low 
mass galaxies akin to COLD GASS, and resolved gas 
maps of galaxies spanning a broader range of properties. 

The main difference between our global results and 



those of studies focussing on the (sub-)kpc scale may 
rather be more methodol o gy-ba sed. As pointed out cor- 
rectly by IRahman et al.l (|2012t ). resolved studies mea- 
sure i<2e P (H2) in regions of the galaxies where CO is de- 
tected, but studies like ours include the total gas mass 
and the total star formation rate. The resolved measure- 
ments therefore teach us about star formation in the gas- 
rich GMCs and the global measurements about the total 
timescale for star formation, including any diffuse gas 
or star formation component. In particular, the global 
values of id e p(H2) include star formation coming from 
the outer disks, regions where CO is in general not de- 
tected in resolved maps and therefore not included in 
resolved studies. Since star formation is known to be 
very inefficient in these outer regions (id e?) (HI)> t Hubble, 
IBigiel et al.l 120101 ) , the relative contributions of outer 
disk star formation to the total value may be respon- 
sible for some of the galaxy-to-galaxy star formation ef- 
ficiency variations, and the difference in absolute value 
of the mean depletion time reported in different studies 
(in general a factor of ~ 2 higher in resolved studies). 
For example, iBolatto et all (|2011| ) find < d ep(H 2 )~ 1.6 
Gyr at the kpc scale in the SMC, but tdep{^2)= 0.6 
Gyr for the entire galaxy, in line with the values we 
find in galaxies with low stellar mass surface densities. 
There are however arguments against the "outer disk ef- 
fect" interpretation. First, there is evidence that the 
molecular star formation law extends even in the HI- 
domi nated outer regions of normal star- fo rming discs 
(e.g. iHever et al.l 120041: ISchruba et all 120111 ). Secondly, 
the metallicity gradients of galaxies in the stellar mass 
range M* > 10 10 M(^ are rema rkably flat out to the optical 
diameter (|Moran et al.l 120121 ) . such that there is no evi- 
dence that the CO measurements of COLD GASS galax- 
ies are subject to radially-dependent aco variations. 

High resolution mapping of CO in a broad range of 
galaxies will be necessary to further make the connec- 
tion between local and global measurements. As the lat- 
ter are most commonly available for high redshift sys- 
tems, understanding this relation is critical. Surveys like 
COLD GASS, offering homogeneously measured quan- 
tities over large, representative samples of galaxies are 
essential to identify the key populations where star for- 
mation efficiency is modulated, allowing us to understand 
the physical processes that are responsible for regulating 
the star formation budget of the local universe. For ex- 
ample, while bars can lead to nuclear starbursts and the 
formation of bulges, they are not responsible for increas- 
ing significantly the star formation efficiency over the 
local galaxy population. Galaxy interactions and merg- 
ers however do have that effect. Now that statistically- 
large samples are sta rting to be collected at high redshifts 
(jTacconi e t al. 2010], and in prep.), it will also be possi- 
ble to track over cosmic time the relative contributions 
of different mechanisms such as bar instabilities, bulge 
stabilization, interactions/mergers and AGN feedback in 
the regulation of star formation. 
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APPENDIX 

A. MORPHOLOGICAL CLASSIFICATIONS 

One of the objectives of this work is to evaluate the impact of interactions on the star formation efficiency, requiring 
us to first identify the subset of COLD GASS galaxies undergoing such events. We do not wish to restrict the analysis 
to systems that can be classified as major mergers or that have massive companions, of the type that could be identified 
in the SDSS spectroscopic database. Instead, we choose to identify all the galaxies showing signs of morphological 
d isturbances. 

I Wang et al.l (|2010l ) measured for all the galaxi es in the GASS/C OLD GASS survey two different quantitative mor- 
phology parameters, based on the approach of lLotz et al.l ([2004). The asymmetry parameter (Asym) provides a 
measure of the difference between the SDSS g— band image and a copy of the same image rotated by 180°. The larger 
the value of Asym, the more asymmetric the galaxy. The second is the dumpiness parameter (Clump), computed 
from the difference between the original g— band image and a copy of the same image smoothed to a scale of 20% 
of the Petrosian radius. Galaxies with clumpier morphologies therefore have larger values of Clump. Two additional 
parameters are commonly used to quantify the morphologies of galaxies: the Gini (G) and M20 parameters. The Gini 
coefficient resembles the concentration index, but does not require the use of aper ture photometry or th e identification 
of the galaxy center, since it instead relies on the distribution of pixel values (Abraha m et al.l [2003). It therefore 
performs well in very clumpy systems with ill-defined centers of light. The M20 index is also akin to con centration, 
but a gain considers the full pixel distribution and therefore does not rely on the use of regular apertures ([Lotz et al.l 
12001 . The SDSS .g-band images are used to measure G and M20 for all COLD GASS galaxies. 
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Figure 8. Different diagnostic plots for morphologically disturbed g alaxies based on the four parameters Clump, Asym, M20 and the Gini 
index. In panel (b), the green line is the cut proposed bv lLotz et al.l IMPOST ) to identify major mergers. Plotted are all galaxies detected in 
CO. The colored symbols represent the galaxies visually identified as being interacting/merging. The red symbols are the galaxies disturbed 
morphologically (with or without close companions), and the orange symbols represent the major mergers. 



To assess the performance of these four indicators at isolating the morphologically-disturbed systems we wish to 
study, we also perform a visual classification. Galaxies are marked as "disturbed" if they satisfy one of three criteria: 
(1) they have a companion within 1' (i.e. 30-60 kpc in our redshift interval of 0.025 < z < 0.050), (2) tidal features 
are visible, or (3) the stellar disks are lopsided or show strong asymmetrical features. The classification was done 
"blindly" , without any knowledge of the value of the morphological parameters or any other physical properties of the 
galaxies that were inspected. AS performed the classification twice, randomizing the order of the galaxies each time 
to avoid "fatigue biases" . Only galaxies consistently classified as disturbed in the two runs are flagged as such. 

In Figure [8] the COLD GASS galaxies with a detection of the CO line are shown in different diagnostic plots that 
combi ne the four morphological parameters (Asym, Clump, G and Af2o)- The M20 — G plane was used bv lLotz et al.l 
(2008) to identify major mergers. Figure|5jD shows that their criterion (G > — O.I4M20 + 0.33) does indeed pick up the 
most extreme systems, but not the full range of morphologies we wish to track. 

The visually-identified disturbed galaxies however cleanly separate from the rest of the sample in the Asym-G and 
Asym-Clump planes (Figs. [HK,c). The Asym parameter appears to be the most sensitive tracer of lopsidedness and 
low surface brightness tidal features in the outer disks. A threshold of Asym > 0.35 was for example suggested by 
IConselice et al.l (|2003f ) to identify merging galaxies. There are a number of outliers that are not visually identified as 
disturbed, but located in Figure [5] above the various selection thresholds. These galaxies tend to have the smallest 
angular sizes (r^o^ < 5") or to have high inclinations (i > 60°), two situations where morphological classifications, 
both automated and visual, tend to be less reliable. 

The analysis presented in §3.31 is based on the visual morphological classification. To confirm that our observation 
of shorter molecular gas depletion times in morphologically disturbed galaxies is not simply the product of our clas- 
sification scheme, we repeat the analysis shown in Figure [5] on the full sample with different selection criteria. The 
results are shown in Figure O We compare the initial result (panel a) with the cases when either the morphological 
parameters are used to select disturbed galaxies (panel b), or when we reject from the analysis the most edge-on 
galaxies, where classifications are uncertain (panel c). In all cases, a KS test confirms with very high significance that 
tdepiJAz) is shorter in disturbed galaxies compared to the rest of the gas-rich population. Even when restricting to 
the subsample with low stellar mass surface densities, no matter how we define morphologically disturbed galaxies or 
whether or not we include edge-on galaxies in the analysis, we always find a KS probability < 5% that disturbed and 
un-disturbed galaxies have similar idej^Eb) distributions. In other words, we confirm that low stellar mass surface 
densities are linked with lower than average £dep(H2), and that within this population, disturbed/interacting galaxies 
have the shortest depletion times. 
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Figure 9. Molecular gas depletion time distributions for various subsamples among the COLD GASS galaxies with a secure detection 
of the CO line (S/N > 5). In each panel, the light gray histogram shows the parent population considered in this specific case, the blue 
solid line the galaxies within this parent population having logfi* < 8.7MQkpc — 2 , and the filled dark green the sub-population classified 
as disturbed/interacting. The specific selection criteria in the three cases are as follows, (a) Disturbed galaxies are identified based on 
the visual classification, (b) Disturbed galaxies are selected in the Asym — Clump plane (Fig. |8ji), (c) The parent sample is restricted 
to galaxies with inclinations < 60°, to avoid the most inclined galaxies where morphological classifications are more uncertain, and the 
disturbed galaxies are identified based on the visual classification. In each case, the KS probability that samples are drawn from the same 
parent population are given for two sets of samples: KS a ;; compares the parent sample and the disturbed subsample (i.e. the gray and 
green filled histograms, respectively). No matter the selection criteria, the depletion times of the disturbed subsample is always significantly 
different compare to the entire gas-rich population. The second value, KSjo^,,, compares the disturbed subsample among the population 
below the critical value of logfi* = 8.7Af0kpc — 2 (i.e. the greed filled histogram and the open blue histogram). 

B. THE CHOICE OF a C o AND ITS IMPACT ON THIS STUDY 

The total molecular gas mass is commonly derived from the measured CO line luminosity using a conversion factor 
olco as Mh 2 = a coL' CO M y where ^co(i-o) * s * n ^ m s P c2 )- The exact value of aco and its possible dependence 
on a range o f parameters including metallicity and density is being active l y investigated, bo t h with observations (e.g. 
Ilsraell Il997t IDame et all 120011: IRosolowsky et ail 120031: [Blitz et alj|200l [Leroy et all 120111: IGenzel et all 12012ft and 
simulations (e.g. iGlover fc Mac Lowll2011t Ishettv et al1l2011t iNaravanan et alJl2012tlFeldmann et al.ll2012ft . 

There is a g eneral consensus that within the Milky Way aco= 4.35 M W (K km s _1 pc 2 ) -1 , with little spatial 
variations (e.g. IStrong fc Mattoxl 119961 : IDame et all 120011 : lAbdo et"aLl|20 iol). There is also evidence for at least two 
sets of circumstances where aco departs significantly from the Galactic value: low metallicity environments and the 
very dense central re gions of merging s ystems . 

We investigated in iSaintonge et all (|2011bft the possibility that depletion time variations observed in the COLD 
GASS sample are caused by metallicity-induced changes in aco- While a few individual galaxies with M*~ 1O 1O M0 
may have a higher-than-galactic aco-, we rule out the possibility that the observed t^ep^i) variations across the 
sample are caused by such aco variations. The various prescriptions for the metallicit y dependence of ag o indicate 
that departu r es fro m the Galactic value become measurable at Z < 0.5Z Q (see e.g. ILerov etHTI ([201 lft . but also 
IGenzel et all ([2012ft f or a stronger ^-depend ence, albeit at z > 1), while most COLD GASS galaxie s have Z > O.8Z 
on the 03N2 scale of IPettini fc Pagell ([2004) . adopting the solar abundance of lAsplund et all (|2004ft . 

For a subset of 196 COLD GA SS galaxies, high qu ality nebular abundances are available from a companion optical 
spectroscopy observing program ([Moran et al.ll2012ft . In Figure [TUl the depletion time distributions are presented in 
three intervals of both luminosity- weighted, slit- integrated metallicity and extinction. The idep(H2) distributions for 
the three metallicity intervals are statistic ally equivalent, as confirmed by a KS test. This is to be compared with for 
example Figure 10 of ISchruba etall (|20T1 : in our stellar mass range (where metallicities are high and their gradients 
flat), the global values of idep(H2) do not depend on metallicity 

Models however suggest that an even better tracer of aco is extinction, as this quantity is the convolution of the 
metallicity and the volume density ([Glover fc Mac Lowll201lft . We find that the galaxies with the highest extinctions 
(Ay > 1-5) have longer tdepi)^2) than the rest of the sample (KS pro b = 0.04 that the distributions are drawn from the 
same parent population), however since these galaxies also have larger than average /i*, this may only be the reflection 
of the previously reported dependence of id e p(H 2 ) on /i*. More interestingly, there is no sign for the least extincted 
galaxies to have shorter than average idep(H2); galaxies in the lowest quartile of Ay have i<2 e p(H 2 ) values consistent 
with the rest of the population (KS pro ;,=0.56) 

Additionally, as shown in Figure [TTT the break in the i^o-SFR (or Leo ~ Lfib.) relation occurs at M*~ 1O 9 ' 5 M , 
below the mass range of the COLD GASS sample. To extend the sample, in particular in the mass range of 
M*< 10 10 M^, which is not prob ed by COLD GASS, Figure [TT] also pr esents data from the lite r ature compila- 
tion o f iGracia-Carpio et all (1201 1|). which includes CO measurements fr o m [Sanders fc Mirabel (11985ft: iMirabe 



([1990ft; ISanders et all (11991ft: ICombes et all (Il994ft: I Young et all (Il995ft: lElfhag et all (11996ft: IMaiolino et all 
Heifer et al l (12003ft: iGao fc Solomon! ti2004ft: lEvans et all (12005ft: I Albrecht et all (12004 12007ft: i Kuno et all ^ 
Lerov et all ((2009) ; Yo ung et all (|2011ftT rG arcia-Burillo et all ([2012ft . A s we further do not obs erve any dependence of 



the residuals of the idep(H2)-M* relation on metallicity (see Fig. Bl oflSalntong e et aI1l2011bft . we conclude that the 
COLD GASS galaxies, bar a few possible exceptions, are not in the regime where aco is affected by photodissociation 



18 



Saintonge et al. 




log t dep (H 2 ) log t dep (H 2 ) 



Figure 10. Left: Molecular gas depletion time in three different metallicity intervals: 12 + log(0/H) < 8.65 (open red histogram), 
8.65 < 12 + log(0/H) < 8.75 (dashed blue histogram) and 12 + log(0/H) > 8.75 (filled green histogram). A KS test indicates that the three 
distributions are consistent with being drawn form the same parent sample, arguing that metallicity effects on aco ar e n ot at the root 
of the observed depletion time variations. Right: Same as left panel, but for three bins of extinction, Ay, as measured from the Balmer 
decrement. The metallicity and extinction bins are defined such as to isolate the top and bottom quartile in each parameter, respectively. 
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Figure 11. Relation between CO(l-O) luminosity and star formation for the COLD GASS sample (open squares) and a compilation of 
galaxies from the literature (filled circles). Symbols are color-coded according to their stellar mass. The break in the relation between 
Lqo ar >d SFR, highlighted by the light-blue shaded region, occurs at M*~ 10 9 ' 5 Mq, below the mass range of the COLD GASS survey 
(Af*> 10 10 M Q ). The histogram shows the SFR distribution of the COLD GASS galaxies with a reliably measured CO(l-O) flux. 

of the CO molecules due to reduced dust-shielding. 

However, some of the galaxies in the sample are interacting, and it has been shown that in syste ms undergoing major 
mergers such as nearby ULIRGs and high— z SMG s, th e value of aco is reduced by a factor of 5 (jSolomon et alJll997t 
iDownes fc Solomonlll998t[Tacconi et al J I2008) . In 32. 31 we have presented our favored strategy to assign a value of aco 
to each galaxy. In short, a galaxy is given either the Galactic or the "ULIRG" aco based on a quantitative measure 
of the ISM conditions, in this case the dust temperature as parametrized by the Sqo/Sioo ratio. As a reference, the 
KS relation derived with this conversion factor choice and initially presented in Figure |5] is repeated in Figure 112b . 
In panel (a) of the same figure, we show the KS relation, had we assigned a Galactic aco for all the galaxies. The 
main result of this study, namely that interacting/disturbed galaxies have shorter than average depletion times while 
bulge-dominated galaxies have the longest depletion times, is still visible. However, in this case the galaxies undergoing 
the most violent interactions (and also having the highest IR luminosities and dust temperatures) do no longer mark 
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Figure 12. The star formation (Kennicutt-Schmidt) relation from the COLD GASS sample for different aco predictions: (a) a Galactic 
value of ctco = 4.35 A/q(K km s — 1 pc 2 ) - 1 is applied to all galaxies, (b) galaxies are assigned either a Galactic aco or the "ULIRG" value 
of 1.0 Mq(K km s — 1 pc 2 ) -1 as explained in £12.31 This is the approach used i n this paper. Finally, in pa nel (c) we use a prescription for 
aco that depends smoothly on Td ua t as parametrized by the Sao/Sioo ration ( Gracia-Carpio ct al. 2011). 

the upper envelope of the distribution. The second consequence of applying the Galactic aco to the merging galaxies 
is that their gas fractions become very large (0.4 < Mh%/M* < 1.0), unrealistically so considering that normal 
star-forming galaxies in our sample have on average Mhi/M* = 0.09, and never more that Mh%/M* = 0.3. 

In Figure [12b. we show the results of applying a pres cription for ago that also depends on Sso/Sioo, but this time 
in a continuous manner rather than as a bimodal one ()Gracia-Carpio et al.ll20lil) . The results are both qualitatively 
and quantitatively similar. 



